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We study the experimental signatures of top flavored dark matter (top FDM) in direct detection 
searches and at the LHC. We show that for a dark matter mass above 200 GeV, top FDM can be 
consistent with current bounds from direct detection experiments and relic abundance constraints. 

We also show that next generation direct detection experiments will be able to exclude the entire 
perturbative parameter region for top FDM. For regions of parameter space where the flavor partners 
of top FDM are not readily produced, the LHC signatures of top FDM are similar to those of other 
models previously studied in the literature. For the case when the flavor partners are produced at 
the LHC, we study their impact on a search based on transverse mass variables and find that they 
diminish the signal significance. However, when the DM flavor partners are split in mass by less 
than 120-130 GeV, the LHC phenomenology becomes very distinctive through the appearance of 
displaced vertices. We also propose a strategy by which all parameters of the underlying model can 
be experimentally determined when the flavor partners can be observed. 


I. INTRODUCTION 

Astrophysical observations have revealed conclusively 
that the dominant type of matter in the universe is a 
particle which is not part of the Standard Model (SM), 
referred to as Dark Matter (DM). While the gravitational 
interactions of DM have been firmly established from an 
observational point of view, several classes of experimen¬ 
tal searches have so far not succeeded in providing any 
information about the properties of DM as a particle. 
In the most commonly considered classes of models, DM 
arises as a thermal relic with weak-scale cross sections 
and the preferred DM mass is not too far removed from 
the masses of SM particles, typically in the GeV-TeV 
range. 

The majority of models for beyond the Standard Model 
physics in this energy range aim to address the Higgs 
mass hierarchy problem by pairing SM particles with 
partner particles of the same gauge quantum numbers, 
but not necessarily of the same spin. Since DM must not 
carry electromagnetic or color charge, there are only a 
few choices in these models for who the DM particle can 
be partners with, and the phenomenologically most vi¬ 
able identities for DM are the partners of the electroweak 
gauge bosons, or neutral degrees of freedom in the Higgs 
sector. 

However, one must not forget that there is an implicit 
assumption in all these models, namely that DM and the 
solution to the Higgs mass hierarchy problem are inti¬ 
mately connected. There is no direct evidence for this 
assumption to be true, which brings up the question of 
whether there are any other guiding principles one might 
adopt in constructing models of DM. If one drops the 
assumption that DM is tied to the solution of the hierar¬ 
chy problem, other interesting approaches present them¬ 


selves, such as a possible connection between all matter 
particles, including DM as well as the SM matter degrees 
of freedom. 

All SM matter degrees of freedom have a property in 
common: they come in three generations. The origin of 
the SM flavor structure is unknown, and it is assumed to 
arise from physics at a high energy scale. If one pursues 
the alternative guiding principle of looking for a connec¬ 
tion between DM and matter particles in the SM, one 
is then led to the possibility that DM itself may trans¬ 
form nontrivially under flavor symmetries. It is therefore 
worthwhile to consider classes of models with “Flavored 
Dark Matter” (FDM) [THIS], to look for universal fea¬ 
tures of such models, and to study whether there exist 
experimental signatures that may be used to distinguish 
FDM from more traditional models of DM, such as a 
neutralino in a supersymmetric setup. It is even possi¬ 
ble that this line of inquiry may yield new insights in 
searching for the origin of the SM flavor structure. 

The general consequences of this approach and the 
classes of models that arise from this starting point were 
considered in Ref. [5j. In particular, the FDM models 
of phenomenological interest contain a local interaction 
that couples the DM flavor multiplet to a SM matter field 
via a mediator, which in the most minimal setup is a fla¬ 
vor singlet. We will refer to this as the FDM interaction, 
and we will use <f> to denote the mediator field. The FDM 
interaction has several consequences that determine the 
general phenomenological features of this class of models. 
First of all, it is natural for the DM flavors to be split in 
mass, because even if the DM flavors have the same mass 
at a high energy scale, mass splittings will be induced at 
loop level from the FDM interaction. Second, the heav¬ 
ier DM flavors can decay to the lightest one, therefore 
only the lightest member of the DM flavor multiplet is 
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cosmologically stablcj^] From this point on we will refer 
to the stable component of the DM multiplet as x, and 
to the heavier flavors collectively as x'■ Third, when <f> 
is a flavor singlet, constraints from flavor violating SM 
processes strongly prefer the DM and SM flavor struc¬ 
tures to be aligned, which means that the DM multiplet 
should have three flavors. These then become flavor part¬ 
ners with the SM matter field that the FDM interaction 
couples them to. The flavor alignment is most easily for¬ 
mulated in the language of Minimal Flavor Violation, or 
MFV [HI, although a more general scenario has also been 
studied I Hi] . 

At this point one can still choose SU( 2) x U( 1) quan¬ 
tum numbers for the DM multiplet as well as which SM 
matter multiplet it is coupled to in the FDM interaction. 
In the interest of minimality, we will restrict ourselves 
to the case where the DM multiplet is a SM gauge sin¬ 
glet. The SM multiplet appearing in the FDM interac¬ 
tion and the specific flavor in that multiplet that couples 
to the lightest FDM state x is then used to label the 
FDM model in consideration, for example in a model of 
electron-FDM, \ is directly coupled to the electron in the 
FDM interaction. There are two possible choices even at 
this point, since the FDM interaction may contain the 
right-handed or the left-handed lepton field. Again in the 
interest of minimality we will choose <j> to be an SU( 2) 
singlet, which means that only SU( 2) singlet SM fields 
appear in the FDM interaction, so that electron FDM 
would couple to the right-handed electron. 

The scenario of lepton-flavored DM has been studied 
in the literature mmm and shown to have many in¬ 
teresting features, including the possibility of having an 
asymmetry between the DM and its antiparticle, gener¬ 
ated during high-scale leptogenesis. For quark flavored 
DM, direct detection bounds can be severe if x is coupled 
to the light quarks, because tree-level (j> exchange leads 
to a large cross section for x scattering off of nuclei [5]. 
On the other hand, when x couples to a third genera¬ 
tion quark, the dominant process for x-nucleus scattering 
arises at loop level, and direct detection bounds are less 
constraining. In this paper, we will turn our attention 
to the case of top FDM, and we will study the current 
exclusion bounds as well as future discovery prospects in 
direct detection experiments as well as in LHC searches. 
We will show that relic abundance considerations com¬ 
bined with direct detection constraints can generically 
be satisfied for m x > 200 GeV. In this parameter range, 


1 Unless the lifetime of the heavier flavors is much longer than the 
age of the universe, which is not generic and requires a rather 
special setup. In this paper we will only concern ourselves with 
the generic case. 


not only do the existing LHC searches at 7 and 8 TeV 
have little exclusion power, but we will show that even at 
the 14 TeV LHC, searches based on traditional kinematic 
observables only become sensitive to top FDM at high lu¬ 
minosity > 100 fb” 1 . When x ' is heavy and is not readily 
produced on shell at the LHC, the collider phenomenol¬ 
ogy becomes similar to other models that have been pre¬ 
viously studied in the literature, in particular Ref. [13] . 
On the other hand, we will show that when the \ and 
X'' masses are split by less than 120-130 GeV, top FDM 
models give rise to displaced vertices that would not only 
make discovery significantly easier, but also allow for all 
parameters of the underlying model to be experimentally 
determined. We will also show that the next generation 
of direct detection experiments will be able to probe the 
entire perturbative parameter region that is not currently 
excluded. 

The layout of the paper is as follows: In Sec. [TI] we will 
review the theoretical aspects of the top FDM model in 
detail. In Sec. [In] we will calculate the cross section for x 
to scatter off of nuclei, and we will determine the region 
of parameter space that has been excluded by existing 
direct detection constraints, as well as the region that 
can be probed with the next generation of direct detec¬ 
tion experiments. In Sec. |IV| we will discuss the LHC 
signatures of top FDM, we will evaluate the prospects 
for a search strategy based on kinematic observables and 
we will study under which conditions signal events can 
be expected to contain displaced vertices. We will also 
comment on how the parameters of the model can be 
measured post-discovery. We will conclude in Sec. [V] 

II. THE MODEL 

In this section we will briefly review the theoretical fea¬ 
tures of top FDM. We will use the conventions of Ref. [5], 
in particular, we will take the DM multiplet to transform 
in the (anti-)fundamental representation of I/(3) u c. The 
Lagrangian contains the FDM interaction 

C D A ijXiUj 0 + h.c. (1) 

Here, \i denotes the DM flavor triplet and <j> is the flavor 
singlet mediator. As mentioned in the Introduction, the 
DM multiplet will be taken to be a SM gauge singlet, 
which means that (f> carries both SU( 3) color as well as 
hypercharge. Note that there is a global 17(1) symmetry 
under which <fi has charge +1, all three have charge —1 
and all SM fields are uncharged. This symmetry can be 
named “x-number” and it ensures that the lightest state 
that has nonzero x-number is stable. Since <j> carries color 
and electric charges, it is not phenomenologically viable 
for it to be stable, so in this paper we will only consider 
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the case that <f> is heavier than the lightest of the \i- 

From the interaction term of Eq. [l] it is clear that one 
of (j> and x must be a fermion while the other is a boson. 
As in Ref. [5], we will adopt a benchmark for the rest of 
this paper that x is a Dirac fermion while <f> is a complex 
scalar. Different setups are possible, such as making \ 
a boson and 4> a fermion, but these do not dramatically 
change the qualitative features of the phenomenology. 

The simplest way to be consistent with flavor con¬ 
straints is to adopt the MFV paradigm [Tf||, namely 
that the SM Yukawa couplings are the only source of 
the breaking of flavor symmetries. Clearly, this is not 
intended to be a complete model. Rather it is an as¬ 
sumption about the UV completion of the SM in which 
the underlying flavor symmetry of the full model is spon¬ 
taneously broken by a single source, from which both the 
SM Yukawa couplings and the FDM coupling of Eq. [I] in¬ 
herit their flavor structure. While it is a very interesting 
problem to construct such a UV model, this is beyond 
the modest goals of this paper, and we will simply adopt 
the MFV paradigm as an assumption here. 

With the MFV paradigm, and the assumption that 
U(3) x is identified with U( 3) u c, the A^ matrix has a fla¬ 
vor structure corresponding to 3 x 3 of U (3) u c. Assuming 
further that the relative normalization of the spurions Y u 
and Yd is preserved by the UV theory, we can then write 

Atf = (at+ 0 YjY u ... + 'YYjVcKMY d Yj V^ KM Y U ...) , 

( 2 ) 


where Y u =diag {y u ,y c ,y t ) and Y d =diag (y d ,y s ,y b ). Note 
that in addition to the first two terms, which are the 
unique contributions to order Y° and Y 2 * with the cor¬ 
rect flavor structure, we have also included the term with 
coefficient 7 , which is only one of several possible terms of 
order Y 4 that has the correct flavor structure. The reason 
we have singled out this particular term is that it is the 
leading term that cannot be diagonalized due to quark 
mixing in the SM. Since y 2 y 2 , for most purposes this 
term can be neglected against the term proportional to 
0, but we will return to it in Sec. |IV C] where naively even 
a small flavor-changing coupling can make a difference, 
if only to show that it can be safely neglected even then. 

Since y u>c <C 3 /*, it is in principle possible for the cou¬ 
pling of Xt to the SM to be different than those of Xu,c 
(which can be taken as equal to each other to a very 
good approximation) by an order-one factor, if 0 ~ a. 
The results of this paper will be essentially independent 
of the fact that these two values may be different, and 
therefore for simplicity we will consider the case 0 -C a 
for most of this paper, in which case all three couplings 
have the same value, which we will denote simply by A. 
In Sec. IV D we will return to this point and consider 


how one might attempt to separately measure the values 
of A t and A UjC at the LHC. 

MFV also dictates the flavor structure of the x mass 
term to have the form 

i m x\ij = ( m ol + A to YJY u ).. . (3) 

Unlike the case of the coupling matrix A ij, here the pres¬ 
ence of the leading correction proportional to A?n has 
great phenomenological significance. As mentioned in 
the Introduction, there cannot be any relic Xu or Xc left 
over today because in that case DM can scatter off of 
nuclei at tree level, and this possibility is strongly con¬ 
strained by direct detection experiments (see Ref. [5j for 
quantitative details). Therefore Xt must be the lightest 
state. Depending on the value of the mass gap, the de¬ 
cay of Xu and %c to Xt may proceed through on-shell or 
off-shell top quarks. We will study the branching ratios 
for this decay in more detail in Sec. |IV| Not surprisingly, 
the value of Am is quite important for collider physics, 
as it will affect the kinematics of signal events containing 
these heavier states. Note also that since y 2 -C y 2 <C y 2 , 
the mass gap between Xu and Xc is much smaller than 
the mass gap between either one of them and Xt, and 
therefore for all purposes of this paper, Xu and Xc can be 
taken to be degenerate. Indeed, when the very slightly 
heavier Xu state is produced, it would decay directly to 
Xt since the available phase space for it to decay to Xc is 
extremely small. As already mentioned in the Introduc¬ 
tion, we will refer to Xt as x from this point on, and both 
Xu and Xc will be denoted by x!■ We will only distinguish 
between them when not doing so can lead to confusion. 

The relic abundance of top FDM is set via the 
usual weakly interacting massive particle (WIMP) mech- 
anis In principle, coannihilations between the three 
X flavors may be important. However, for the reasons 
outlined above, it is preferable to take Am to be appre¬ 
ciable, and we will work with the simpler scenario where 
all x! have already decayed before x freezes out. As will 
be shown in Sec. m direct detection experiments favor 
m x > 200 GeV , therefore the dominant annihilation 
channel is XX _■* U through <fi exchange, and the annihi¬ 
lation cross section (s-wave) is given by 


_ n__ X _ 

ailcl 327t (m x + — mf) 2 

While we will not pursue this possibility further in 
the rest of the paper, it is worth remarking that in the 



2 As proposed in Ref. ESI, the flavor structure of FDM makes it 
possible to generate an asymmetry for \ i n the early universe. 

This intriguing possibility for top FDM is beyond the scope of 

this paper, but will be taken up in a future study. 
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case when % is lighter than the top quark, the freeze- 
out calculation is significantly modified, both because the 
top quark drops out of thermal equilibrium significantly 
above m x , and also because the annihilation has to pro¬ 
ceed through off-shell tops, or in extreme cases, through 
the channel XX > 99 via a top quark loop. 


III. DIRECT DETECTION CONSTRAINTS 


to leading order in m x /m ( j ) bjj^] 
3A 2 e 


9x = 


m 2 t 


327T 2 swcw 
3\ 2 eQ t 


1 + log ^4 
mi 


4 ~ ~ I ^ 3 l0 S to 2 


Mx = 


647r 2 m0 

3A 2 eQ t m x 
647T 2 TOi 


( 6 ) 


The amplitude for top FDM to scatter off of a nucleus 
in a direct detection experiment receives contributions 
from a number of processes, all arising at loop level since 
X is a SM singlet. In particular, at loop level \ develops 
charge-radius and magnetic dipole type couplings to the 
photon and the Z boson (shown in Fig. [l]J^] a coupling 
to the Higgs boson (essentially the same process as the 
diagram on the left of Fig. [T] but with 7 /Z replaced 
by h), and effective four-point couplings to u-u and g- 
g generated through box diagrams. These contributions 
were all calculated in Ref. [5], with the conclusion that 
the 7 /Z exchange diagram dominates over the four-point 
coupling to u-u , the remaining diagrams being negligible 
compared to these two. 



FIG. 1. Feynman diagrams that contribute to the effective 
coupling to 7 /Z in the top FDM model. The vector boson 
lines in the loop diagrams can be attached to either the top 
quark or to the mediator <f> running in the loop. 


Let us start by focusing on the diagrams of Fig. |T| 
After integrating out <j>, the one-loop induced couplings 
of top FDM are given by the effective Lagrangian 

£eff = 9 X XTIgX^ + b x XlvXd^ u + HxXiVuvXF^, (5) 

where the effective coupling g x of the Z boson, the charge 
radius b x and the magnetic dipole moment /i x are given 


Using Eq. [5] we can write down an effective Lagrangian 
describing the scattering of x from a light quark q as 

k 

£eff = c|x7 M X97/x 9 + c®X7 /i X97M<7 + 

(7) 


where 

4 = 4 = e< 3A> 4d = eQqkx’ (8) 

rn z 

and g q z is the coupling of the quark q to the Z boson. 

We next convert the quark-level operators to nucleon- 
level operators and take the non-relativistic limit, arriv¬ 
ing at the following effective Lagrangian at the nucleon 
level (N = p , n) 

£ e ff = cfXb^xAhv-ZV + c^xYx^l^N 
+ c QXi° a ^xNK, J N 

+4xicr^xNia^kpN, (9) 

where K fL is the sum of the incoming and outgoing nu¬ 
cleon four-momenta. Here the coefficients c N are given 
by 


jv _ 9z9x 

C Z — 

m z 

Cry - CQ 


( 10 ) 

( 11 ) 


with g p 7 = — 2 —(i — 2s?^) and g z = — 5 —. The 

VZj SwC-W V 2 VV / JZ 2swC\V 

coefficients Cq and are given by 

Cq = eQ N /j, x /2m N , c™ = -efl N p, x /2m N , (12) 

with flp = 2.8 and jl n = —1.9. 

So far we have kept the magnetic dipole terms. The 
momentum dependence in the magnetic dipole interac¬ 
tion makes it impossible to factor the differential event 
rate as a product of the elastic cross section and the ve¬ 
locity integration. Following the procedure in appendix 


3 We drop magnetic dipole couplings to Z since those are sup¬ 
pressed by fc 2 /m| where k is the typical momentum exchange 4 For all plots in this paper, we use the exact expressions of g x , b x 
in the scattering. and q x which are also valid for m x < m^,. 
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FIG. 2. The LUX bound on the coupling A for = 500 GeV 
calculated using the effective Z coupling and charge term, 
with and without the addition of the dipole term. 


C of Ref. [IS], we calculate the differential rate numer¬ 
ically, and work out the exclusion limits from LUX [2D] 
in the presence of the dipole terms. The result is shown 
in Fig. [2] We find that the effect of the magnetic dipole 
operator is negligible compared to the effective Z and 
charge operators in setting limits for the coupling A. 
Based on this conclusion, for the rest of the paper we will 
neglect the magnetic dipole contributions, which will al¬ 
low us to directly use the LUX bounds for the elastic scat¬ 
tering cross section in order to evaluate the constraints 
on the parameter space of the top FDM model. 

In a similar fashion we study the importance of the 
effective four-point coupling to u-u by using the result of 
the box diagram calculated in Ref. [5]. In Fig. [3] we plot 
our results in the form of the LUX bound on the coupling 
A with and without taking the box diagram into account, 
confirming the conclusions derived in Ref. [B] that the di¬ 
rect detection bounds are dominated by the effective j/Z 
coupling. For the rest of this paper, we will use only this 
dominant contribution in calculating the bounds from di¬ 
rect detection experiments on the parameter space of the 
top FDM model. 

The leading (spin-independent) contribution to the 
matrix elements at the level of nucleons is 

(x, N X, N) = Am x m N . (13) 

It is then convenient to define the dark matter-nucleon 
effective couplings 

C N = 4ni x rriNCz + 4m x mNC^ , (14) 

in terms of which the total scattering cross section (at 



FIG. 3. The LUX bound on the coupling A for = 500 GeV 
calculated using the effective Z coupling and charge term, 
with and without the addition of the effective coupling to 
u-u. 


the level of the nucleus) can be written as 

<7 T = -*-f ---) [ZC P + (A - Z)C n ] 2 . (15) 

167 T \m x + m p J 

We show in Fig. [4] the LUX [2D] exclusion region in the 
A -m x plane (A -m^ plane) for a generic value of (m x ). 
For TO 0 = 500 GeV, A is constrained to be smaller than 
0.7. Note also that the correct relic abundance (shown 
by the blue curve) is obtained for m x > 200 GeV. This is 
generic: Independent of the value of when m x gets 
close to the top mass, the relic abundance constraint fa¬ 
vors large values of A which is in conflict with direct de¬ 
tection bounds. Also, for any fixed value of m x , obtaining 
the correct relic abundance puts an upper limit on m$ as 
can be seen in the second plot of Fig. [4] 

In Fig. [5] we show the LUX exclusion region in the 
m x~ m 4’ pl ane for fixed values of A. Not surprisingly, the 
exclusion region grows with A. However, the curve on 
which the correct relic abundance is obtained also ex¬ 
tends further, such that an order one fraction of the pa¬ 
rameter space still remains viable. In order to illustrate 
this even more efficiently, in Fig. [bj we plot the (green 
shaded) LUX exclusion region by assigning to each point 
in the mass plane the value of A for which the correct 
relic abundance is obtained for that mass point. As can 
be seen from this plot, only relatively light x masses are 
excluded. Contours corresponding to different values of 
A are also plotted for reference. 

Having mapped out the exclusion region from current 
constraints, we turn our attention to what region of pa- 
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m x [GeV] 



mt [GeV] 


FIG. 4. The LUX exclusion region in the A -m x plane (A-m^ 
plane) for m$ = 500 GeV (m x = 200 GeV). The correct relic 
abundance is obtained on the blue curve. 


rameter space can be probed with the next generation of 
direct detection experiments. We use the projected ex¬ 
clusion reach of the XENON1T experiment with an expo¬ 
sure of 2.2 ton years [21] . which is provided in the range 
m x < 1 TeV. We find that top FDM (subject to relic 
abundance constraints) can be entirely excluded in this 
range. This is shown in Fig. [6] as the blue shaded region. 
We also extrapolate the XENON1T projected exclusion 
reach to m x > 1 TeV by using n x oc m~ 1 and assuming 
a constant detector efficiency for m x < 10 TeV. We find 
that the exclusion region extends to the full perturbative 



[GeV] 

FIG. 5. The LUX exclusion bounds in the plane for 

A = 0.5 (green shaded region) and A = 0.6 (blue shaded 
region). For each case, we also show the curve on which the 
correct relic abundance is obtained. 


parameter range (A < 1) for the model. 


IV. DISCOVERY PROSPECTS AT THE LHC 

We now proceed to investigate the discovery prospects 
for top FDM at the LHC. We will focus on the case 
m x > 200 GeV, which as we have shown in the previ¬ 
ous section is strongly favored by direct detection exper¬ 
iments and relic abundance constraints. Searches based 
on kinematic features such as missing transverse energy 
or transverse mass variables have not been very ef¬ 
fective in probing this mass range at the end of the 7 TeV 
and 8 TeV runs of the LHC, and therefore we will focus 
our attention on discovery prospects at the 14 TeV run. 


A. Production and Decay Modes 

Since x is a SM singlet while 0 is colored, the most im¬ 
portant production mechanism is 0-pair production from 
gluons in the initial state. This production mode is pure 
QCD and is therefore independent of A. 0 can also be 
pair produced through QCD interactions from a q-q ini¬ 
tial state. However, in this case there is also a contribu¬ 
tion from t-channel x' (specifically, \u) exchange, which 
has A dependence. Processes contributing to 0-pair pro¬ 
duction are shown in Fig. [7] and the cross section for 
several choices of A are shown in Fig. [8] We see that 
the t-cliannel x' exchange can have an order one effect 







































7 



[GeV] 

FIG. 6 . The current LUX exclusion bounds in the 
plane (green shaded region), where A is varied to give the 
correct relic abundance at each point. Contours of constant 
A are also shown for reference. The full parameter region in 
this plot (and beyond) can be excluded by the XENON1T 
experiment (blue shaded region). 


on the cross section for larger values of A, which are not 
excluded by direct detection bounds. 




FIG. 8 . (f >-pair production cross section for A = 0, 0.5,1 with 
m x ' = 300 GeV. 



FIG. 9. Feynman diagrams for <i>-X associated production 
and y'-pair production at the LHC. 



FIG. 7. Feynman diagrams for the pair production of the 
mediator (f> at the LHC. 


()>-pair production is not the only production channel in 
this model, since y„ can couple to the initial state quarks 
through the FDM interaction of Eq. jT| In particular, 
both 0-Xu associated production and y M -pair production 
are allowed, and these proceed via the Feynman diagrams 
of Fig. [9j Note that cr(pp -0 fiXu ) oc A 2 while a(pp -0 
XuXu) oc A 4 . The cross sections for these channels are 


shown in Fig. 10 


Note that all production channels involve </> and y u , 
both of which decay promptly. Therefore both sides of 
each event will contain visible final states, cj) has “short” 
and “long” decays, where in the former case it decays 
directly to a top quark and y, while in the latter case cj) 
decays to a light quark and the corresponding flavor of y', 
which subsequently decays to a light quark, a top quark 
and a y through an off-shell q i. Depending on the mass 
difference between y' and y, the top quark itself may be 
off-shell. This is illustrated in Fig. mi The branching 
ratio for the short decay mode is shown in Fig. |T2] as a 
function of the ratio m x /m x ' for several values of rn x '. 
The branching ratio for the long decay mode (which is 
really two decay modes for y„ and y c , which have iden¬ 
tical branching ratios to a very good approximation) can 
simply be found as one minus the branching ratio for 
the short decay mode. The production and decay modes 
under consideration ultimately result in a t-t pair and a 















4)X’ Production Cross Section, A=1 





<t> V & x. 

FIG. 11. The short and long decay modes for cf>. 


X'x' Production Cross Section, A=1 



m^’ (GeV) 


FIG. 10. (Top) 4>~x' associated production cross section for 
A = 1 and two values of m x >. a scales like A 2 . (Bottom) 
y'-pair production cross section for A = 1 and two values of 
m$. a scales like A 4 . 


X~X pair i n every event (the tops possibly being off-shell 
depending on the details of the spectrum), and anywhere 
between zero and four additional light jets depending on 
the number of on-shell x! hr the event. 

Note that while m x < must be satisfied in or¬ 
der to obtain a neutral DM particle, if Am in Eq. [3] is 
large enough, then m x > can exceed m^. When that is the 
case, the collider phenomenology is particularly simple: 
all signal events basically arise from 0-pair production 
followed by the short decay mode of 0. The LHC discov¬ 
ery prospects of this simpler case are essentially the same 
as for the model studied in detail in Ref. m • Therefore 
in the remainder of this section we will focus our at¬ 
tention on how the collider prospects are affected by the 
presence of y 7 , and what kind of impact this may have on 
the conclusions reached in Ref. 13]. Note that the answer 
is not entirely obvious. On the one hand, one may argue 



FIG. 12. Branching ratio of the short decay mode of 0 as a 
function of the ratio m x /m x > for several values of m x ’ . When 
m$ — m x drops below rtit, the short decay mode can only 
proceed through off-shell tops and therefore the branching 
ratio becomes negligible. 

that the presence of x' may improve discovery prospects 
by opening new production modes, thus increasing the 
cross section, and by making the long-decay mode of 0 
possible, which gives rise to additional light jets in the 
event. On the other hand, for values of A that are allowed 
by direct detection bounds, the collider phenomenology 
is still dominated by 0-pair production, thus the number 
of signal events does not significantly increase, and fur¬ 
thermore within these events, those with the now acces¬ 
sible long 0 decays may be more difficult to distinguish 
from background, since the top quarks and the $ T in 
these events are only produced in the last decay stage, 
which due to the presence of x! has a smaller available 
phase space, therefore leading to softer tops and less $ T . 
Finally, the presence of x! may allow for a more exotic 
signature in signal events, namely displaced vertices that 
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occur in x! decays, which may be more important than 
both of the above considerations for discovery. This last 
possibility however may or may not be realized depend¬ 
ing on the mass spectrum. 

In the following, we will first attempt to evaluate how 
the presence of x' affects the discovery prospects for top 
FDM in a search based on transverse mass variables. We 
will do this by setting up a relatively simple analysis 
based on the dileptonic decay mode of the top quarks, 
which is the cleanest channel. We will look at several 
mass spectra for a fixed value of with varying m x 
and rn x ', since kinematic features of signal events will 
depend on the differences of m^, m x > and m x , but they 
will be less sensitive to the overall mass scale. In partic¬ 
ular, as the mass differences in the spectrum change, top 
quarks produced in decays may be on or off shell, leading 
to very different efficiencies for analysis cuts. In contrast, 
the total signal normalization will depend mainly on the 
overall mass scale (m^), but it will be insensitive to the 
mass differences. 

After studying the impact of x! i n a kinematics-based 
search, we will return to the question of displaced ver¬ 
tices in x' decays, and then comment on how these may 
allow all the parameters of the top FDM model to be 
determined experimentally. 


B. Search Based on Transverse Mass Variables 

Event Generation and Selection: Since we will fo¬ 
cus on the dileptonic decays of the top quarks, the signal 
events in our analysis are characterized by two 6 -jets, 
a pair of opposite sign leptons that are uncorrelated in 
flavor, and $ T . SM processes with similar final state 
particles constitute backgrounds for the top FDM signal. 
In order to model the relevant backgrounds, we gener¬ 
ate samples of 500,000 events (2 million for t-t) at 14 
TeV center of mass energy for the backgrounds listed 
in Table [T] The t-t background samples are generated 
using MC@NLO (22, 23j whereas the remaining back¬ 
grounds are generated using MadGraph [23] at LO, and 
later rescaled with A'-factors. The order of the A-factors, 
and the associated references are also listed in Table |U 

We force all final state tops and W’s to decay lepton- 
ically, including leptonic decays of r’s. The ttZ back¬ 
ground has a much smaller cross section compared to 
the others. However, we include it because our search 
strategy will be based on Mt2 [28H5D] and ttZ it is the 
dominant background at high Mj- 2 - For this reason we 
force the Z boson to decay to neutrinos, as any other 
decay mode does not contribute at high Mt 2 and is sub¬ 
dominant at low Mt 2 - 

In generating signal samples, we fix m$ = 500 GeV. 


The production cross section plots of the previous sub¬ 
section can be used to extrapolate the results of the anal¬ 
ysis to different values of The mass points we use are 
given in the first two columns of Table [TT] For each value 
of m x we also include a spectrum with x' decoupled as 
a point of comparison in order to better assess the im¬ 
pact of the presence of on discovery prospects. This 
is not unrealistic, since even for my moderately above 
TO 0 , 0 -pair production will be the dominant production 
mode and essentially no on-shell x! will be produced. In 
order to be consistent with relic abundance and direct 
detection constraints, we only consider m x > 200 GeV, 
and we fix A = 1/2. Note that this choice of A is conser¬ 
vative for the 0-pair production cross section (see Fig.| 8 |, 
and decay kinematics are insensitive to A. For each mass 
spectrum we generate 200,000 events at 14 TeV using 
MadGraph [23] for each of the signal production modes, 
which we combine with the appropriate ratios. All top 
quarks in signal samples are forced to decay leptonically. 

Both signal and background samples are passed 
through FastJet El E2 to reconstruct jets using the 
anti-A:* jet clustering algorithm with radius R = 0.4 [33] . 
We use a 6 -tagging efficiency of 0.6. Isolated leptons are 
identified with cuts on the variable pt( iso) around the 
lepton track, where the pr( iso) is the sum over the pt 
of the tracks with pt > 1 GeV for each track, within 
a A R < 0.2 cone around the lepton candidate. For 
electrons, the isolation cut is chosen to be pr(iso) < 
0.1 x p-r(lepton), while for muons, the isolation cut is 
Pt( iso) < 1.8 GeV. Any jet overlapping with an isolated 
lepton within A R < 0.2 is eliminated from the jet candi¬ 
dates, and electron or muon candidates overlapping with 
any remaining jets within A R < 0.4 are discarded to 
avoid accepting leptons from the decay of heavy hadrons. 

Jet energy resolution is simulated using a Gaussian 
smearing with a variance a: 


a 

E 


S N 


(16) 


where the terms S,N,C are chosen to be 0.5 GeV 1 / 2 , 

5 GeV, and 0.04, respectively. The rj and 0 of the four- 
momentum of a jet is also smeared by a Gaussian with 
a = 0.025. Electron energies are smeared with a = 0.02 x 
Y^A/TeV and muon energies by a = 0.1 x \JE/ TeV. 

We also smear the transverse missing energy of the $ T , 
defined as the magnitude of the two-vector momentum 
summed over all the final state invisible particles (neu¬ 
trinos and LSPs). The individual x and y components of 
the ftjrp vector are smeared by [34] 

a($ T ) x , y [GeV] = (0.40+0.09 x0S)x\/ H T [GeV] + p x 20, 

(17) 

where Ht is the scalar-sumnred pt of all visible final state 
particles within ( 77 ! < 5.0, and p denotes the number 
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Background 

Order 

<7 after cuts (pb) Reference 

tt 

NLO 

14.3 

m 

tw 

NNLO 

0.7 

m 

tiZ 

NLO 

2.3 x 10~ 3 

m 


TABLE I. Details of background simulation. The t-t background is generated at NLO, while the others are generated at LO 
and rescaled with A-factors. The order of each A'-factor we use and the associated references are listed, as well as the cross 
section of each background after selection cuts. 


of interactions per bunch crossing. For each event p is 
chosen from a Gaussian distribution with central value 
equal to 50 and standard deviation equal to 10. 

For the analysis, we use the following selection cuts in 
order to reduce backgrounds: 

• At least two jets with pt > 20 GeV and \rj\ < 2.5 

• At least one ^-tagged jet 

• Two opposite-sign leptons (e and p only) with p t > 

20 GeV and |? 7 | < 2.5. Events with same flavor 
dileptons are vetoed if | mu — mz \ < 10 GeV 

• > 20 GeV. 

The cross section of each background sample after selec¬ 
tion cuts is listed in Table [Tj 

Analysis and Results: Our search strategy will rely 
on the kinematic variable Mt 2 [25H5D] in the leptons + 
]firp part of the event, which is theoretically bounded from 
above by rriw hr events where the missing transverse en¬ 
ergy arises entirely from two neutrinos originating from 
W decays, whereas for events with additional sources 
of missing transverse energy, Mt 2 can be larger, and is 
therefore a good variable to separate signal from back¬ 
ground. Mt 2 is defined as[^] 

max 2 )) 

1 rT,l ' r t, 2 rT J 

(18) 

where the minimization is performed over all possible 
partitions of the missing transverse energy into two parts, 
subject to p T + 9 = p . This models the neutrino 

transverse momenta in the event. mr(p T ,£i) is the 
transverse mass of the subsystem consisting of the hy¬ 
pothetical neutrino and the respective lepton, which for 
massless particles is defined as 


Mj 1 2 — 


mm 

-4-ri \ 


After the selection cuts described in the previous sub¬ 
section have been applied, we look for the choice of an 
Mt 2 cut for each mass spectrum that will maximize 
S/y/B, keeping the number of signal events after all cuts 
above 50 at 100 fb -1 of integrated luminosity. In Ta¬ 
ble [XT] we list the spectra we use, the cross sections for 
each production mode after selection cuts, and the opti¬ 
mal Mt 2 cut for each spectrum, along with the resulting 
value of S/y/B. A plot of the Mt 2 spectrum for the 
background and a generic signal point (m x > = 400 GeV, 
m x = 200 GeV) is shown in Fig. 

These results are not difficult to interpret. The power 
of the Mt 2 based analysis arises from the production of 
X particles with appreciable momenta. It is easy to see 
from Table |Tl] that the most important factor for discov¬ 
ery prospects is the mass gap between <j> and y, as this 
sets the upper limit of the \ momenta, in particular in the 
short decays of <fi which yield the most energetic x among 
all production and decay channels. We also see that the 
value of m x / is the second-most important quantity that 
impacts discovery prospects. The best cases in terms of 
S/y/B are obtained with the spectra in which x/ is de¬ 
coupled. This shows that the significance of the search 
is driven by the short decay mode of (f> which produces 
more energetic tops and more $ T . When x' is reintro¬ 
duced with a mass less than m^, it affects the results of 
the analysis in two ways: 1) For fixed and m x , as m x > 
increases, the branching ratio for the short decay mode 
increases, which boosts S/y/B. 2) Once m x >-m x goes 
above the top threshold, on-shell tops can be produced 
in the long decay mode of <j> as well as in decays of / 
in the <j)-y/ and x'-x' production modes. This boosts the 
signal significance as well. 
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C. Displaced Vertices 


”i!(Pt,i,Pt, 2 ) = 2 Ipt.i||Pt,2|(1 - cos Av?p Til , PT2 ). 

(19) 


5 In our analysis we only use A/7'2 with zero input mass for the 
invisible particle. 


The lesson from the search based on transverse mass 
variables appears to be that the presence of / hurts 
rather than helps discovery prospects for top FDM. Thus 
the discovery reach mapped out in Ref. |13j . which cor¬ 
responds to top FDM with / decoupled, is the best case 
scenario, while for > m x ' > m x , the events with 
on-shell / have less JS T and less energetic tops. What 
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m x i (GeV) 

m x (GeV) 

07 , 4 , (fb) 

(fb) 

a xx (fb) Best M!$ 

(GeV) S/VB 

N/A 

300 

7.7 

- 

- 


< 1 

400 

300 

9.0 

3.4 

0.11 


< 1 

N/A 

250 

8.1 

- 

- 

145 

3.3 

400 

250 

9.0 

4.8 

0.29 

135 

2.0 

350 

250 

9.0 

4.1 

0.14 

130 

1.7 

N/A 

200 

8.4 

- 

- 

165 

8.9 

400 

200 

9.0 

4.8 

0.36 

155 

5.7 

350 

200 

9.0 

5.3 

0.39 

145 

3.5 

300 

200 

9.0 

4.9 

0.19 

145 

3.5 


TABLE II. The signal spectra used in our analysis, with = 500 GeV. For each spectrum we list the cross section after 
selection cuts for each production mode, with A = 0.5. We also list the choice of the Mti cut that maximizes S/VB, subject 
to the constraint that at least 50 signal events pass all cuts at 100 fb -1 , as well as the value of S/\/~B for that choice. For 
the spectra with m x = 300 GeV, S/xB remains below 1 for any choice of Mt 2 cut. For each value of m x we also include a 
spectrum where x' is decoupled as a point of comparison in order to better assess the impact of the presence of on discovery 
prospects. 


I T2 spectrum of SM backgrounds and Top FDM signal 

)6 A —.—.—.—i , , ,—.—i—.—, , , r —■—■—■—i , , ,—■—i—■—, , , 



FIG. 13. Mt 2 histogram of background and the signal point 
(m$,m x i ,m x ) = (500,400,200) GeV with A = 0.5. The ver¬ 
tical line shows the best Mt 2 cut. Background histograms 
are stacked while the signal spectrum is shown separately. 


on-shell x! do add to signal events is light jets. However, 
it is difficult to find an optimal search strategy relying on 
light jets, both because light jets are also easily produced 
in the background from additional initial and final state 
emissions, and also because the pt spectrum of light jets 
in the signal will be highly dependent on the mass dif¬ 
ferences in the top FDM model and it is therefore not 
possible to devise cuts that will be optimal over the full 
parameter region. 

There is however one case where on-shell y; 7 , either di¬ 
rectly produced via the <j)-x! and x'~x' production modes, 
or arising from the long decay mode of </>, may signifi¬ 


cantly improve discovery prospects. When m x ’ — m x < 
m t , the x' decay becomes four-body (x! bWjx)- The 
lifetime of x! is then in the range for displaced vertices. 
This is shown in Fig. 14 If m x ' — m x is even lower, 


namely below mw +mb, then even the W boson must go 
off shell, and the x! decay becomes five-body. While in 
the four-body regime the decay length cr can be resolved 
with the silicon vertex detectors, in the five-body regime 
cr becomes truly macroscopic. For at least a part of the 
five-body regime, one would expect to observe x' decays 
occur in the calorimeters or even the muon system. 


These conclusions assume that the FDM interaction 
of Eq. |T] is diagonal in flavor space. However, we know 
that even if the UV theory respects MFV perfectly, there 
will be flavor-changing decay channels that y 7 will inherit 
from quark mixing in the SM. This is in fact why we 
chose to include the term proportional to 7 in Eq. [2j as 
it allows us to estimate the size of the leading contribu¬ 
tion to flavor-violating x! decays in the exact MFV limit. 
Fortunately, even for 7 ~ 0(1), the flavor-violating (tree- 
level) three-body decay Xc —l► cc\ will be suppressed by 
(ylycV c bY which makes it subdominant to the flavor¬ 
preserving decay even in the five-body regime. A fla¬ 
vor violating decay of \u would of course be even more 
suppressed, and therefore we conclude that the MFV- 
respecting corrections in Eq. [2] do not affect our con¬ 
clusions about displaced vertices. Loop-induced flavor- 
violating decays are even smaller, for instance Xc —■► Xtl 
with cf> and a charm quark running in the loop can be ob¬ 
tained from the above described tree-level flavor-violating 
decay by contracting the external charm quark lines and 
attaching a photon to the loop. The coupling structure 
in flavor space remains the same and therefore leads to 
the same suppression. Moreover, even though the decay 
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is now two-body, the added suppression from the loop 
factor and a™ ~ 0( 10~ 2 ) makes the partial width for 
this decay channel even smaller than the tree-level flavor- 
violating decay. 

In the absence of a complete UV model, we should of 
course only think of MFV as an approximation. However, 
deviations from MFV cannot be of order one in order 
to avoid tension with bounds from SM flavor-changing 
processes. Taking into account the suppression from the 
smallness of any off-diagonal 0 couplings dictated by such 
SM constraints, the predictions for the four-body decay 
mode should be unaffected. On the other hand, whether 
or not the five-body decay mode can ever be dominant 
will depend on exactly how small such beyond-MFV off- 
diagonal couplings are. Therefore, when discussing the 
possibility of macroscopically displaced decays, we 
should not forget that there is a significant amount of 
model dependence involved. 

Note also that there are loop-induced processes that 
allow x' decay without any flavor violation even for 
very small m x > — m x . The leading such process is shown 
in Fig. [15] At first sight, this process may appear to be 
of the same size as the tree-level decay in the five-body 
regior^ However, a closer inspection reveals this to be 
incorrect: both quark propagators in the loop are at¬ 
tached to an FDM vertex on one end and to a IF boson 
on the other. However, while the FDM interaction in¬ 
volves purely right-handed quarks, the weak interactions 
of course couple purely to left-handed quarks. Therefore 
both quark propagators in the loop require a mass inser¬ 
tion. This is not a big price to pay for the top quark, 
however already for the charm quark, this mass inser¬ 
tion ensures that the loop induced decay remains sub¬ 
dominant to the tree-level decay, even in the five-body 
regime. 

D. Measuring the Parameters of the Model 

In the rest of this section we will comment on how the 
parameters of the top FDM model can be measured once 
discovery level significance has been achieved. As our 
conclusions so far have made clear, the answer to this 
question will be very different in various regions of pa¬ 
rameter space of the model. Below, we will go through 
the distinct possibilities, starting from the case where cer¬ 
tain parameters simply cannot be determined, towards 
more optimistic cases where all model parameters can be 


Based on the very rough rule-of-thumb that the suppression from 
adding a loop is equivalent to the suppression from adding two 
final state particles at tree level. 



FIG. 14. The decay length cr of x' as a function of m x i for 
fixed values of = 500 GeV, m x = 200 GeV and A = 0.5. 
The results at very low mass splitting should be regarded as 
an idealization only, since in this regime deviations from MFV 
can cause flavor-changing x' decays to become the dominant 
decay channel. The decay length changes by many orders of 
magnitude as m x i — m x goes through the mt and mw + mb 
thresholds. The range in which the LHC experiments are 
expected to be most sensitive to displaced decays is indicated 
as a gray-shaded horizontal band. 



FIG. 15. The leading contribution to x! decay which is flavor 
preserving, without a top quark in the final state. 

measured. Note that while the procedures we will out¬ 
line below can work in principle, they will likely require 
a high-statistics signal sample, and whether or not this 
condition can be met will depend on the details of the 
model, in particular the 0 mass which is the primary fac¬ 
tor in setting the overall production cross section. The 
full list of model parameters is the mediator mass m^, 
the DM multiplet masses m x and m x / (or equivalently 
mo and Am of Eq. [ 3 ]), and the couplings A t and A„ iC 
which may in principle be different (or equivalently a 
and (3 of Eq. [ 2 ]). 

Case 1) x' decoupled: As discussed earlier, if m x ' > 
m 0 , essentially no x' are produced on shell. The LHC 
phenomenology is dominated by 0-pair production and 

































































































13 


furthermore for the values of A* that are allowed by di¬ 
rect detection and relic abundance constraints, the pro¬ 
duction in this channel is dominated by QCD processes 
(see Fig. [8]). This means that experimentally m x ' and 
A u ,c play no role anywhere, and therefore cannot be mea¬ 
sured. 

Since we are considering measurements that will be 
performed once discovery level significance has been 
achieved, the kinematic information present in the data 
can be used for measuring the mass difference — rn x . 
such as the $ T spectrum or the endpoint of the Mti dis¬ 
tribution of the event including the 6-jets (rather than 
the Mt 2 observable we used in the analysis presented 
earlier in this section which was constructed from the 
leptons and ]fi T alone). Combined with the signal pro¬ 
duction cross section, which is a function of alone, the 
two masses TO0 and m x can then be independently de¬ 
termined. The parameter X t does not affect the collider 
phenomenology. However, once the <f> and % masses have 
been measured, one can make a prediction for the value 
of At based on the relic abundance constraint. Further¬ 
more, in the event that a signal is observed in a direct 
detection experiment, the predicted value of At can be 
cross-checked with the signal event rate in that experi¬ 
ment. 

Case 2) m x ' < m 0, m X ' — m x > m* : Signal events 
in this case will sometimes contain on-shell y', either di¬ 
rectly produced, or originating from the long decays of 
cj>. However, the mass splitting between x! and \ is large 
enough that x! decays are three-body and prompt. The 
measurements discussed in case 1 can still be performcc0 

In order to go any further, one needs a robust way 
to separate out signal events with and without on-shell 
x'■ Since there are no displaced vertices, and since top 
quarks produced in x' decays are on-shell, there is no 
simple way of doing so based on kinematics, although 
one may rely on the existence of additional light jets in 
the event which occur in events containing at least one 
x'■ If requiring additional hard jets allows one to extract 
a signal subsample containing on-shell then a detailed 
kinematic study of x' decays, combined with a knowledge 
of 77i0 and m x determined as in case 1 can then be used 
to infer the x! mass. 

As long as x! decays are prompt, kinematics alone will 
not yield any information on the final remaining param¬ 
eter A MjC . However, this coupling does have phenomeno¬ 
logical importance: The production cross section for the 


cj)-x' channel scales as A,„ c , the cross section for the xl~xl 
channel scales as X^ c , and the short vs. long branching 
ratios of 0 depend on X UtC as welQ One can therefore 
attempt to measure X UtC by separating the signal into 
bins with different numbers of light jets, and performing 
a global fit to the number of signal events in each bin us¬ 
ing the predicted production cross sections and branching 
ratios in all channels as a function of X U:C . 

Case 3 ) m x i < m 0, m x / — rn x < mt but no dis¬ 
placed vertices: Note that the x! decays do not be¬ 
come displaced immediately as the x'~X mass splitting 
drops below the top quark mass. In fact, depending on 
the value of A, the x' decay length only becomes large 
enough to be detected by the LHC experiments when 
the mass splitting is 120-130 GeV (see Fig. 141. How¬ 
ever, even though displaced vertices are not observable, 
there is a method based on kinematics in this case that 
allows one to obtain a subsample with a higher fraction 
of events containing on-shell x' ■ This relies on the fol¬ 
lowing observation originally made in Ref. |35j : the px 
spectrum of 6-jets produced in on-shell top quark decays 
peaks at a value that is independent of the boost of the 
original top quark. When the top quark is off-shell how¬ 
ever, the 6-jet pt spectrum is peaked towards zero. This 
suggests that by imposing a suitably chosen upper cut on 
6-jet pt’s, one can reject more events with no on-shell x ' 


than those with on-shell x! ■ This is illustrated in Fig. 16 


After obtaining a y'-enriclied subsample of events this 
way, one can then attempt to measure m x > and A„ iC in a 
similar way as described for case 2. 

Case 4 ) m x > < m 0, with displaced vertices: When 
the mass splitting between x! and X is less than 120- 
130 GeV, the x' decays contain displaced vertices. This 
makes it much easier to obtain a signal subsample with 
on-shell y', and determine m x > by studying the kinemat¬ 
ics of x! decays. Furthermore, one now has a much more 
robust handle on A„ :C , since it can be measured directly 
from the decay length of x! ■ 


V. CONCLUSIONS 

We have studied the phenomenology of top FDM at 
direct detection experiments and at the LHC. In order 
to obtain the correct relic abundance, the preferred mass 
of the DM particle should be 200 GeV or higher. This is 
entirely compatible with current bounds from direct de¬ 
tection experiments. We have also studied the prospects 


7 Note that even though several production channels may be con¬ 
tributing to the signal, the events with the highest Mx 2 still come 
from the short decay mode of </>, so the Mx2 endpoint yields the 
same information as in case 1. 


In principle even the production cross section in the 0 -</> channel 
has a nontrivial A u ,c dependence, however if A u ,c < 0 . 5 , this 
additional contribution is subdominant and therefore difficult to 
observe. 
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FIG. 16. The enrichment factor, defined as the enhance¬ 
ment of the fraction of signal events containing at least one 
x' after demanding that both 6-tagged jets in the event have 
Pt < X, with X plotted on the x-axis. The red curve 
is obtained from a signal spectrum with my — m x < rtit 
while the green curve is obtained from a signal spectrum with 
m x >—m x > mt. As described in the main text, the x' content 
of the former signal point can be enhanced using such a cut 
with a low X because x' decays with off-shell tops give rise to 
softer 6-jets, whereas for the latter signal point with on-shell 
tops, the b-jets tend to be harder. 

of next generation direct detection experiments to look 
for top FDM and we have shown that the full pertur¬ 
bative parameter region (A < 1) will be probed by the 
XEN0N1T experiment with 2.2 ton years of exposure. 

There are several production and decay channels for 
top FDM at the LHC which populate final states with 
two (possibly off-shell) top quarks, missing transverse en¬ 
ergy and possibly additional light jets. In the case that 
the heavier DM flavors x' are not readily produced at the 
LHC, the collider phenomenology is simple, and similar 
to that of other models that have been studied previ¬ 
ously in the literature such as the model of Ref. m- 
We have therefore focused mainly on the differences with 
this simpler case when x! does play a role in the phe¬ 
nomenology. In particular, we have shown that when the 
DM flavors are sufficiently split in mass, the presence of 
X 1 actually reduces the discovery reach of searches based 
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on transverse mass variables. However, when the mass 
splitting is less than 120-130 GeV, we have shown that 
the decay length of x! will be in a range that leads to 
displaced vertices. In this case, signal events are much 
easier to separate from background and the presence of 
x' significantly improves discovery prospects. 

We have also commented on how the parameters of 
the underlying model can be experimentally determined. 
Even in the case when x! are not produced on shell, a 
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tion and relic abundance constraints should allow one 
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